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Abstract. We present self-consistent dynamical models for dust driven 
winds of carbon-rich AGB stars. The models are based on the cou- 
pled system of frequency-dependent radiation hydrodynamics and time- 
dependent dust formation. We investigate in detail how the wind prop- 
erties of the models are influenced by the micro-physical properties of 
the dust grains that enter as parameters. The models are now at a level 
where it is necessary to be quantitatively consistent when choosing the 
dust properties that enters as input into the models. At our current 
level of sophistication the choice of dust parameters is significant for the 
derived outflow velocity, the degree of condensation and the estimated 
mass loss rates of the models. In the transition between models with and 
without mass-loss the choice of micro-physical parameters turns out to 
be very significant for whether a particular set of stellar parameters will 
give rise to a dust-driven mass loss or not. 



1. Introduction 



Mass loss by dust driven winds of asymptotic giant branch (AGB) stars prob- 
ably is one of the major mechanism which recycle material in the Galaxy (e.g. 
Sedlmayr 1994). Most stars {M^ < 8Mq) will eventually become AGB stars and 
subsequently end their life as white dwarfs surrounded by planetary nebulae. 

AGB stars are cool (T* < 3500 K) and luminous (L* of a few 10^ to a few 
10^ Lq) and a majority of them are pulsating long-period variables (LPVs). The 
outer layers of many AGB stars provide favorable conditions for the formation 
of molecules and dust grains. Dust grains play an important role for the heavy 
mass loss of these stars (up to (M ~ lO~^M0/yr). 

Pulsation causes an extended atmosphere where the dust is condensing. The 
dust absorbs the light of the central star and re-radiates it at longer wavelengths 
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Figure 1. A schematic drawing of the atmosphere of an AGB star and 
the most important physical processes. (Courtesy of Christer Sandin). 



in the infrared rage (A > 2/xm). The density of material in the circumstellar 
envelope may be so large that the star becomes completely obscured in the 
optical rage. Due to observational difficulties in determining the mass loss rates 
of these objects considerable effort has been put into a theoretical description 
of the mass loss of AGB stars (e.g. Fleischer et al. 1992; Hofner & Dorfi 1997). 
The overall goal is to develop a mass loss description that can be used as input 
to evolutionary models and thereby probe the chemical evolution of the Galaxy. 

We present here how the wind properties of carbon-rich models are influ- 
enced by the micro-physical properties of the dust grains. 



2. The Atmosphere Model 

The models predict mass loss rates, wind velocities and dust-to-gas ratios by 
solving the equations of continuity, motion and energy for the gas and a set 
of equations describing the formation of dust grains. The transfer of radiation 
is calculated frequency-dependently for the gas and dust (Hofner 1999; Hofner 
et al. 2002). Basically only four stellar parameters are needed as input for the 
models; the stellar mass (-M+), the effective temperature (Tefj), the luminosity 
(L*) of the star and the gas abundance ratio (ec/^o) of carbon to oxygen. 
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However, to describe the pulsations of the star two additional parameters are 
required for the prescription of the inner boundary condition; a pulsation period 
(P) and a velocity amplitude (Af), which simulates a pulsating stellar surface 
as a means of mechanical energy input. 

The number densities of the molecules relevant to the carbon dust formation 
are calculated assuming chemical equilibrium between H, H2, C, C2, C2H and 
C2H2 after the fraction of carbon bound in CO has been subtracted. Nucleation, 
growth and destruction of dust grains are supposed to proceed by reactions 
involving C, C2, C2H and C2H2. 

Dust is formed by a series of chemical reactions in which atoms or molecules 
from the gas phase combine to clusters of increasing size. The molecular com- 
position of the gas phase determines which atoms and molecules are available 
for the cluster formation and grain growth. This can lead to a layered structure 
of the circumstellar envelope. Dust formation alters the atmospheric structure 
which significant infiuence the local and global circumstellar envelope dynamics 
and the detailed time-dependent spectral appearance. 

In the case of pulsating atmospheres of Mira and other LPVs the conden- 
sation of dust grains is strongly connected with the shock wave structure of the 
circumstellar envelope, usually triggered by a 2-step mechanism. First there is 
nucleation of critical clusters followed by growth to macroscopic dust grains. 
Grain growth will usually proceed far from equilibrium and it is necessary to 
use a time-dependent description of grain growth to determine the degree of 
condensation and other relevant properties of the dust. 

The grain growth is treated by the moment method which describe the time 
evolution of an ensemble of macroscopic dust grains of various sizes and requires 
the nucleation rate as input (Gail &: Sedlmayr 1988; Ganger et al. 1990). The 
moment method needs the following dust properties as external input: (1) the 
optical properties of the relevant dust types, (2) the intrinsic density of the dust 
material(s), (3) the sticking coefficient which describes the efficiency of growing 
dust grains from the atomic precursors, and (4) the smface tension of the diist 
grains which controls the efficiency of nucleation. The purpose of this work is to 
investigate in detail how the predicted wind properties of the carbon-rich AGB 
models are influenced by the choice of these quantities. 

For the mass loss models it is important to know how the uncertainty in the 
chosen dust parameters affects the obtained results. The model parameters are 
listed in Tables 2 — 5. For a given set of stellar parameters and elemental abun- 
dance, the mass loss rate, the outflow velocity and the degree of condensation 
is obtained by averaging the outflow of the time-dependent models. Previous 
models presented in earlier papers (e.g. Andersen et al. 1999) were all based 
on gray radiative transfer while the models presented here are calculated using 
frequency dependent radiative transfer for the gas and dust (see Hofner et al. 
2002 for details). Dynamic atmosphere models are undergoing a transition from 
qualitative modeling of physical processes to realistic quantitative predictions 
and interpretations of observations. 
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Figure 2. Wavelength dependence of the dust extinction efficiency 
(Qext/fl), see Table 1 for dust annotation. Also shown is a constant 
value (gray opacity) normalised to the Rouleau data at l//m. 



Table 1. List of the different dust data used. 



Reference 


Material 


Pdust 




Designation 




name 


(g/cm^) 


% 


in this paper 


Jiiger et al. 1998 


cel400 


1.435 


67 


Jager 400 


Jager et al. 1998 


cellOOO 


1.988 


80 


Jager 1000 


Rouleau & Martin 1991 


AC2 


1.85 




Rouleau 
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3. Amorphous carbon dust 

Amorphous carbon particles are considered to be a very good candidate for the 
most common type of dust present in circumstcllar envelopes of carbon-rich 
AGB stars. The infrared data of late-type stars generally show a dust emissivity 
law ((5(A) ~ A~^) with a spectral index of /? ~ 1 (e.g. Campbell et al. 1976; 
Sopka et al. 1985; Martin k Rogers 1987; Giirtler et al. 1996). A A'^ behavior 
can be expected in a very disordered two-dimensional material like amorphous 
material (e.g. Huffmann 1988; Jager et al. 1998). Graphite formation in AGB 
stars is ruled out, because of the absence in the observed spectra of the narrow 
band at 11.52yum and the overall shape of infrared graphite spectra which are 
proportional to A~^ (e.g. Draine & Lee 1984). 

Amorphous carbon covers by nature a broad span of material properties. 
The difference in optical properties for three different types of amorphous carbon 
can be seen in Fig. 1. The dust properties needed as input for the models are 
not always know from laboratory experiments. In the case of amorphous carbon 
it has therefore been custom to use values of the intrinsic density for graphite 
material instead of the real value. If the intrinsic density of graphite is used 
instead of the correct laboratory value for amorphous carbon it will influence 
the efficiency of dust formation and therefore the predicted mass loss and colors 
as seen from Table 3. 

4. Results and Discussion 

A frequency-dependent treatment of the radiative transfer in the dynamical 
models has turned out to be crucial for obtaining realistic structures, synthetic 
spectra, and reliable mass loss rates (Hofner et al. 2002). How the choice of 
the micro-physical parameters can change the predictions for the mean outflow 
velocity of the gas and dust, the degree of dust condensation and the mass loss 
is discussed in the following focusing on the optical dust properties, the intrinsic 
dust density, the sticking coefficient and the surface tension. 

4.1. Dust opacity 

The direct influence of the extinction efficiency (Qcxt) on the structure and the 
wind of the dynamical models are shown in Table 2. The models are dependent 
on the choice of laboratory measurements of the extinction efficiency (Qext) as 
already show for gray models in Andersen et al. (1999). In contrast to the 
gray models then for the frequency-dependent models both the absolute value 
and the slope of the dust opacity data become relevant. The absolute value 
of the extinction efficiency has a profound influence on the determined outflow 
velocity of the models, while the slope of the extinction efficiency determines 
the grain temperature. In the models where the Jdger 400 data is used for the 
extinction efficiency of the dust, no dust is formed, since the slope of the material 
opacity dictates a high grain temperature which make the grains unstable against 
evaporation. 
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Table 2. Influence of the extinction efficiency of the dust. Model 
parameters: = 1.0 Mq, sticking coefficients oq = 0.37 = 0.34 
0(C2H = 0.34 ac2H2 = 0.34, abundances are solar except for C/0 = 1.4, 
period P = 650 d, piston velocity Aup = 4km/s. Input: Luminosity 
L^, temperature T^, dust opacity data ^dust) intrinsic dust density pdust- 
Output: Mass loss rate M, mean velocity at the outer boundary (n), 
mean degree of condensation at the outer boundary (/c). 



Model 








Pdust 


M 


(u) 


(/c) 




[Le] 


[K] 




[g/cm^] 


[Mo/yr] 


[km/s] 




113dj 10/9199 


13000 


2700 


Jdger 1000 


1.99 


5.59 • lO"'' 


14.8 


0.05 


113drou/jl85 


13000 


2700 


Rouleau 


1.85 


4.87 • 10-6 


7.38 


0.10 


113dj04pl44 


13000 


2700 


Jdger 400 


1.44 








110djl0pl99 


10000 


2600 


Jdger 1000 


1.99 


7.0 • 10-6 


16.3 


0.10 


110droupl85 


10000 


2600 


Rouleau 


1.85 


2.26 • 10-6 


3.64 


0.12 


110dj04pl44 


10000 


2600 


Jdger 400 


1.44 
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Figure 3. The colors {H — K) and {K — L) for the models (+) pre- 
sented in Tables 2 — 5 and for carbon Mira with a comparable mass 
loss rate; V Hya (A), R Vol (*), R For (x), R Lep (o). 
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4.2. Intrinsic dust density 

We have calculated models with different intrinsic densities for the given material 
assumed to be formed in the stellar wind. When the moment method was 
developed by Gail 8z Sedlmayr (1988) the intrinsic density for the amorphous 
carbon used in that generation of models was not know. The value of graphite 
(2.25 g/cm^) was therefore assumed. This value was later used in almost all of 
the existing models (e.g. Fleischer et al. 1992; Hofner 8z Dorfi 1997). The models 
become much redder when the intrinsic density of graphite is used instead of 
the value for amorphous carbon since significantly more dust is formed in these 
models. Even a small increase of about 10% in the intrinsic density of the dust 
material (as is the case from Model 113djl0pl99 to 113djl0p225) results in a 
doubling of the degree of condensation and a substantial increase in the outflow 
velocity, (u). 

Table 3. Influence of the intrinsic density of the dust. Model param- 
eters: Same as Table 2. 



Model 








Pdust 


M 


{u) 


(/c) 






[K] 




g/cvo? 


[Mo/yr] 


[km/s] 




113djl0pl99 


13000 


2700 


Jdger 1000 


1.99 


5.59 • 10-*^ 


14.8 


0.05 


113djl0p225 


13000 


2700 


JdgerlOOO 


2.25 


7.26 • 10'^ 


20.8 


0.11 


113drou/9l85 


13000 


2700 


Rouleau 


1.85 


4.87 • 10^6 


7.38 


0.10 


113drou/9225 


13000 


2700 


Rouleau 


2.25 


8.24 • 10"^ 


18.1 


0.31 


113dj04pl44 


13000 


2700 


Jdger 4-00 


1.44 








113dj04p225 


13000 


2700 


Jdger 400 


2.25 


2.10 • 10-^ 


1.38 


0.13 



The influence on the mass loss rate is very sensitive to the discrepancy 
between the value of the intrinsic density of the amorphous carbon and the 
value for graphite. For the models 113droupl85 and 113droup225 where the only 
difference is a 20% difference in the value used for the intrinsic density, the 
estimated mass loss rate differ almost by a factor of 2. For the models using the 
Jdger 4OO material data (113dj04pl44 to 113dj04/9225), where the difference is 
almost 40%, using the true material value instead of the higher value for graphite 
result in a model that will not develop a wind at all. 

4.3. Sticking coefficient 

The sticking coefficient (also called the reaction efficiency factor) a is not def- 
initely know as long as we do not know explicitly the sequence of chemical 
reactions responsible for the dust formation. To isolate the importance of how 
uncertainties in the value of a will influence the results of the models we have 
varied this parameter in otherwise identical models. Choosing a value for the 
sticking coefficient of one increases the degree of condensation and the outflow 
velocity by about a factor of 2 compared to when a value of 0.3 is chosen. 
Dust formation occurs beyond « 2R^ independent of the value assumed for the 
sticking coefficient since the condensation threshold is determined mainly by 
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the temperature. With the low choice of sticking coefficient (0.3) only about 
10% of the condensable carbon material present in the gas actually condense 
into grains, while for the much more favorable choice of sticking coefficient (1) 
twice as much material condense into grains. A complete condensation of carbon 
grains is prevented by the rapid velocity increase after the onset of avalanche 
nucleation and the subsequent rapid dilution. 

Table 4. Influence of the sticking coefficient for the dust formation. 

Model parameters: Same as Table 2, except for the sticking coefficient 
0(C C(C2 OiC2H (y-CiH'i.i the dust opacity Atdust = Rouleau, and the intrinsic 
dust density pdust = 1-85 g/cm^. 



Model 






ac 




ac2H 


(^C2H2 


M 


(n) 


(/c) 






[K] 










[Mo/yr] 


[km /s] 




113a03 


13000 


2700 


0.37 


0.34 


0.34 


0.34 


4.87 • 10-*^ 


7.38 


0.10 


113a02 


13000 


2700 


0.20 


0.20 


0.20 


0.20 


3.38 • 10-6 


3.87 


0.09 


113a05 


13000 


2700 


0.50 


0.50 


0.50 


0.50 


5.75 • 10-6 


10.5 


0.12 


113al0 


13000 


2700 


1.00 


1.00 


1.00 


1.00 


7.02 • 10-6 


16.5 


0.22 


110a03 


10000 


2600 


0.37 


0.34 


0.34 


0.34 


2.26 • 10-6 


3.64 


0.12 


llOalO 


10000 


2600 


1.00 


1.00 


1.00 


1.00 


6.95 • 10-6 


12.3 


0.23 



4.4. Surface tension 

The surface tension o"dust of amorphous carbon is not know from laboratory 
experiments, it has therefore become a custom to use the values for graphite 
e.g. given by Tabak et al. (1975). One problem however is that there are huge 
variations in the values determined due to the anisotropy with orientation of 
graphite. 

The surface tension enters into the surface contribution to the free energy, 
AF, calculated on the basis of classical thermodynamics, see Gail et al. (1984) 
for details. 

As already demonstrated by Tabak et al. (1975) then variating the value 
of the surface tension produces an enormous change in the nucleation rate. To 
determine how significant the prescribed value for the surface tension that we 
use in the models is for the results we have varied the value around the value of 
1400 erg/cm^ which were used in all previous models, see Table 5. Altering the 
value of the surface tension of the forming dust grains by 28% around the value 
of graphite will make the difference between obtaining a mass loss or not for the 
model. Also it has a substantial influence on how much of the available material 
in the circumstellar envelope will condense into dust grains. For a comparison 
the measured surface tension for other materials are: Fe (udust = 1400 erg/cm^), 
MgS ((Tdust = 800 erg/cm^) and SiO (cJdust = 500 erg/cm^) (Gail & Sedlmayr 
1986). 
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Table 5. Influence of the surface tension for the dust formation. 
Model parameters: Same as Table 2 but also including the surface ten- 
sion (Jdust aiid with the dust opacity Kdust = Rouleau, and the intrinsic 
dust density pdust = 1-85 g/cm^. 



Model 






CJ'dust 


M 


(n) 


(/c) 






[K] 


[erg/cm^] 


[Mo/yr] 


[km/s] 




113droupl85(7l0 


13000 


2700 


1000 


8.54 • lO-'' 


3.52 


0.76 


113droupl85o-14 


13000 


2700 


1400 


4.87 • 10-6 


7.38 


0.10 


113droupl85(7l8 


13000 


2700 


1800 









5. Conclusions 

We have investigated in detail how the predicted wind properties of the carbon- 
rich AGB models are influenced by the choice of micro-physical dust parameters. 

For the mass loss models it is important to know how the uncertainty in the 
chosen dust parameters affects the obtained results. The choice of the micro- 
physical parameters can change the predictions for the mean outflow velocity of 
the gas and dust by a factor of 4, the predicted degree of dust condensation by a 
factor of 10 and the predicted mass loss by a factor of 2. In the transition between 
models with and without mass-loss the choice of micro-physical parameters is 
vital for whether a particular set of stellar parameters will give rise to a dust- 
driven mass loss or not. 
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